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ABSTRACT 

An axial crack in a cylindrical shell is investigated by use of 
a 10th order shell theory, which accounts for transverse shear defor- 
mations as well as a special kind of orthotropy. The symmetric prob- 
lem is formulated in terms of two coupled singular Integral equations, 
which are solved numerically. The asymptotic membrane and bending 
stress fields ahead of the crack are found to be self similar. Stress 
intensity factors are given as a function of the shell parameter for 
various values of the ratio crack length to shell thickness. Consider 
able differences from 8th order shell theory results are found for the 
bending stresses, while the membrane stresses of the 8th order theory 
seems to be a lower limit reached for very thin shells. 


Part of the present work was carried out during a visit to Lehigh 
University supported by NSF under the Grant ENG 73-045053 A01 and 
NASA-Langley under,, the Grant NGR 309-007-011. 
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1. INTRODUCTION 


Plate and shell theories based on Kirchhoff'.s assumption only en- 
ables the satisfaction of two plate-type boundary conditions. Broadly 
speaking the importance of this restriction depends on the variations 
of the field quantities compared to the plate or shell thickness. 

Thus serious shortcomings can be expected in problems where steep grad- 
ients are encountered, e.g., crack problems. A number of investiga- 
tions [1-4] concerned with a crack in an infinite plate in bending have 
revealed considerable differences between the 4th order bending theory 
solution and solutions obtained by use of 6th order Reissner-type 
bending theory [5], The differences concern both the magnitudes of 
the stresses and their distribution around the crack. In the light 
of these results it seems to be of considerable Interest to supplement 
existing results for cracks in shallow shells described by 8th order 
shell theory with calculations, which explicitly incorporate the effect 
of transverse shear. 

In the present paper a brief derivation of shallow shell field 
equations including transverse shear is given. The procedure follows 
closely that of Naghdi [6], but a certain type of orthotropy is in- 
cluded here. It is demonstrated, how these orthotropic equations can 
be obtained from the isotropic equations merely by use of suitable 
variable transformations. This also holds for the 8th order shell 
theory used in previous investigations [7-10], and these solutions 
can therefore be given a more general interpretation. 



Although the shallow shell equations are fairly general, the rest 
of the paper is confined to an axial crack in a cylindrical shell. 

The method of solution is that of singular integral equations obtained 
by the complex Fourier transform. An investigation of the asymptotic 
stress field around the crack tips is included, and extensive numer- 
ical results are given. 

2. FUWPAMENT AL EQUATIONS 

In this section a brief derivation is given of the fundamental 
equations for a shallow elastic shell. A special kind of orthotropy 
is accounted for through the parameter 5. In the isotropic case the 
equations are those given by Naghdi [6]. The middle surface of the 
shell is described in the cartesian coordinate system of Fig. 1 as 
Z=Z(X a ), a=l,2. The stress resultants are given in terms of their 
components in the {X a ,Z> system. A different approach has been used 
by Si h and Hagendorf [11] for an isotropic spherical shell. 

When only a vertical load q(X a ) is included, the equilibrium equa- 
tions are 


' 0 


V c,a + (Z ,A S >, e + t > = » 


M - V = 0 
n aB,P v a 


( 2 . 1 ) 

( 2 . 2 ) 

(2.3) 


N M n and V are the membrane forces the moments and the transverse 
ap cip _ cl 

shear forces. 


(2.4) 


For a shallow shell the strains are defined by 


where U and W are the displacement components in the {X„,Z> system. 
The normals to the shell in its original configuration change direc- 
tions by the angles 3 . The slope of the middle surface changes by 
W , and thus the effect of the transverse shear is expressed by 

I ^ 


0 


a 



+ 3, 


a 


(2.5) 


From (2.4) a compatibility equation is extracted in the form 

e ccY e 0S^ e a3,Y<5 + Z ,a3 W ,y6^ = 0 
where e is the permutation symbol. 

When elastic shells are considered, Hooke's law yields 


( 2 . 6 ) 


he a0 ' a a8y6 N yfi 


(2.7) 


where h is the thickness of the shell. By use of the stress function 

F(X 1 defined by 
11 


N a3 = e ay e 36 F ,y6 


( 2 . 8 ) 


the equilibrium equation (2.1) is satisfied, while (2.2) takes the 
form 

M a3,aB + Z ..a3 B ay e 3'S F »^ + q = ° * 

The compatibility equation (2.6) becomes 
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(Z.10) 


e aK e 3X e 'Yy e 6v a aeY6 F »KXy\) + hZ ,aB e at e B6 W y5 = 0 


In general solutions of the system of differential equations (2.3), 
(2.9) and (2.10) will be quite complicated. However, considerable simp- 
lifications arise, if the differential operators can be factorized. Be- 
low only orthotropic materials with (X a ) as principal axes are considered. 
This implies the following form of (2.7), 


_ 

e n 


r i 
*7 

V 1 

"*T 

0 


" N ll" 


e 12 

_ 1 
" K 

V 2 

i 

e 2 

0 


N 22 

(2.11) 

2e 12 


0 

0 

1 

Sz_ 


N 12 



where VgE-j = v-|E 2 . In this case the factorization property amounts to 

1 ( 2 . 12 ) 


!i+ 1 
E 1 ^2j 




Now introduce the geometrical mean values 
E = /E^ , v = /v-| \>2 
and the orthotropy parameter 6 defined by 



In terms of these three parameters (2.11) takes the form 


(2.13) 


(2.14) 
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e ll 


&~ z -V 
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” N ll" 

e 12 
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-v <5 2 
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2(l+v) 
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(2.15) 


Following from an assumption of linear variation of the stress compon- 
ents over the thickness the moments are then given by 


H n 


6 2 

v 0 



M 22 
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— — 


A linear relation between the angles 0 a and the shear forces V Q is 
assumed. In order to make elimination of the moments from (2.9) possi- 
ble this relation must be of the form 



(2.17) 


B is the effective transverse shear modulus t which in the isotropic case 
may be taken as I- G [5]. 


The system of differential equations can now be simplified by elim- 
inating the moments. When the notation 


n2 . i 0 ? ' . r- 1 3 Z 

v 6 " 6 axT + 6 m 


(2.18) 
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Is adopted, the result is 


rjl»C 4 . I»P / S 2 Z 8 2 

V t|£| i5fW 


- 2 


a z z r . a'z 3‘ i 
sXISXr 3X753(1 + IXf W' 


(2.19) 


and 


T2" TCF 7 S W " ( ff v o {ffif ' z W^ifafc 


h 2 E , 3 2 Z l! 


3 2 Z 3 s 


+ Mw )f = 


h 2 


(2. 2D) 


with the extra conditions 


&1 + M” * 12(1 -v 2 7 f [7 5 e i + V 3X 2 ( 3Xi ' 3X 2 )] 
h + B TS'CTv 7 ! S’ [V 6 B 2 + 6 *T* ( ^2 (2 ' 22) 


( 2 . 21 ) 


3. PIMENSI0NLESS PARAMETERS 


The following curvature measures are introduced 


1 3 2 Z X - <£?. 

" "3X^ ‘ R z “ " 3X| * R 


]_ 

12 


3 2 Z 


9X 1 3X 2 


(2.23) 


When dimensionless variables and parameters are defined as shown in 
Table 1, the differential equations (2.19)-(2.22) reduce to 


--p- (A 2 - 2X 2 z + A 2 -|jar)w = 0 


(3.2) 


v“w + A 2 (1~kV 2 )(A 2 -g^r “2A 2 2 -^+ A 2 ^t)$ =AM1-<V 2 ) | q 


(3.3) 
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and 


(1-kV 2 )$ x + f£ = K "T 1 ^ (lx” ■ ly 1 ) (3 * 4) 

{l-KV 2 )3 y + |*f « k ^ ^ (*~^ - -§^) (3.5) 

The usual static and geometric boundary conditions can also be 
formulated without explicit use of the orthotropy parameter 6 . As a 
consequence the following presentation refers to an Isotropic shell. 

The corresponding specially orthotropic solutions are easily found 
by use of Table 1. It should be pointed out that the shell parameter 
Xz is 6 Z times the shell parameter used by Yuceoglu and Erdogan [12]. 

4 . THE CYLINDER - INTEGRAL REPRESENTATION 

When considering the cylinder shown in Fig. 2, X i =X I2 =0, The 
homogeneous equations corresponding to (3.2) and (3.3) then take the 
form 


v'* - (^) Z fjx = Cl (4.1) 

V*w + (M^O-cV*) 0= o (4.2) 

Elimination of either w or $ from (4.1) and (4.2) leads to the same 
8th order differential equation 

+ X$(1-k7 2 ) = 0 (4.3) 

Although the introduction of a finite transverse shear stiffness 
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represented by 1/tc does not Increase the order of the differential 
equations (4.1) and (4.2), the order of (3.4) and (3.5) Is Increased 
thereby enabling the satisfaction of five boundary conditions as com- 
pared to four, when k b 0. 

The procedure now Is to give integral representations for $(x,y), 
w(x,y), 0 (x,y) and 3 (x,y). When a crack Is present along part of 

A J 

the x-axIs as shown in Fig. 2, different expressions must be given for 
the half-planes y>0 and y<0. Introduce the representations 


t(x,y) s -Jff f° 5U,y)e' l5x d? 

ioq 

(4.4) 

w(x,y) = ^ £ w(e,y)e“ i?x d? 

(4.5) 

where 5U,y) and w(|,y) are given by 


5(5. y) * ]$(x,y)e 1Cx dx 
*00 

(4.6) 

w(5,y) = f w(x,y)e i?x dx 
^00 

(4.7) 


Due to (4.3) and the corresponding equation for $(x,y) both w(?,y) 
and $(5,y) can be expressed as linear combinations of expOr^y), where 
rn j~ m j (^) are t * ie roots characteristic equation 

(m 2 -5 2 ) 4 + (X^) 4 [l - <(m 2 -e 2 )] = 0 (4.8) 

The notation n . 

p = m 2 - 5 2 (4.9) 
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Is now Introduced leading to the quartlc equation 


p" + (M) 4 {1 -kp) * 0 


(4.10) 


The solutions PjU)» j=1,...4 to (4.10) are given In Appendix A. The 
Solutions nijU) to (4.8) are selected such that 


Re[njj]<0 » m j +4 = " m j » J a ^»-*»4 


(4.11) 


When the displacement transform function wU.y) Is given In the 


form 


■* n^y 

I M5)e J * y>0 

.j=l a 

w(?,y) =<• m.y 

l MS)e 5 , y<0 

U=5 j 


(4.12) 


the transform function i(?,y) Is found by substitution of (4.12) Into 
(4.1). 


4(5»y) =<, 


t \ z 4 r 2 ' m.y 

f- ( f> i, } 


- (?) I (f) R,(S)e 

k X j=5 3 

Introduce the function Y(x,y) as 
^(x.y) = ly- “ *3^ 


V 


y>0 

y<0 


(4.13) 


(4.14) 


The equations (3.4) and (3.5) give the equation 
( 1 - = 0 


(4.15) 


from which the transform function ¥(§,y) is found to be 
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where 




* + 


K\ f-vj 


(4.17) 


Now f(x»y) is assun»ed to be independent of $ x and The equa- 
tions (3.4) and (3.5) are then solved, and the solution verified by 
substitution Into the original equations. When a function x(*»y) with 
the transform 

.V , y> 0 

(4.18) 

» y<o 


x(?»y) ■< 


* EX„ >* 


l i~) 

j*l 3 
# EX, •* l 


V 


Is Introduced, the solution is found to be 


a - + e 1-V. W 

3X + 3y 


8 |I 

p y 5x 2: 3x 


(4.19) 

(4.20) 


The solution to the field equations has now been represented In 
terms of ten unknown functions Aj(5), j*l,2 and Rj(5), j=1...8, which 
must be determined from the boundary conditions of the problem. 

5. B0UMPAW CONDITIONS - SYMMETRIC LOADING 

For the axial crack shown in Fig. 2 the following five static 
quantities are prescribed at the crack surface, 



N yy ' 

■ 14 

ax 7 

• N xy- 

a 2 * 

axay 

K yy ' 

* FT 5 " 

3B 

( — X- + V 

1 ay 

— ) 
ax ' 

M ° 
xy 

If 

l-v , 3fS x 
T ' ay 

30 yx 

av 

Ih., 

3B V 

B i M i 4 

a 2 w 


ax 

ax 

axay 



The last expression is used in differentiated form in order to secure 
dimensional consistency. By substitution of the Integral representa- 
tions of thr, pwviotjs section the following expressions are found for 
the half-plane y>0. The corresponding expressions for y<0 follow trivi- 


ally from changes of the indices. 

if 00 x 2 4 •> m,y -l£x 

Nyy(x.y) ■ ^ [ <r(f) l Pj- Rj(5)e 3 }e de (5.2) 

if” x. z * •> m.y -15x 

« xy (x.y) ■ h [W<T) J,?i R j (5|e J 16 d5 (5 ' 3) 

Hyy(X,y) = j [pj+(l . v) ^ ]pj -4 RjU)e V 

+ igr^U)/ 1 ^ i?X dC (5.4) 

M xy (x,y) = F 7F j 

r y -i£x 

+ 5 J 5 - [1+K(l-v)5 2 ]A,U)e 1 }e d t (5.5) 
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(5.6) 


5 , r * „ ™ ( y 

6 V*-y> * s Vf V )e J 

+ K Iji 5 J A 1 (5)e riy )e' U V 

The problem of static boundary conditions on the crack In Fig. 2 
could In principle be formulated directly in terms of the unknown func- 
tions Aj(S) and RjU). This would lead to a system of coupled dual 
Integral equations, which is not easily solved even by numerical meth- 
ods. An alternative and numerically more suitable technique consists 
In the following. The Integral representations given In Section 4 are 
used to express geometric quantities, which are In a sense complemen- 
tary to the static quantities (5.1). The most direct choice Is the 
generalized displacements corresponding to the generalized forces (5.1). 

For the sake of brevity we shall limit our attention to the case 
of symmetric self-equilibrating loading. 

Nyy(x.y) * N yy (x,-y) » N xy* X,y * s - N xy{*»' y ) 

Myy(x,y) = M yy (x,-y) , M xy (x,y) = -M xy (x,-y) (5.7) 

f ‘ 

V y (x,y) = -V y (x,-y) 

In this case ? 

« i+4 U) = RjU) , a 2 (C) = -A 1 (g) (5.8) 
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whereby the number of unknown functions reduces to five. 


Due to (5.7) the boundary conditions for N , M vu and V are homo- 

xy xy y 

geneous. We therefore only need to Introduce generalized displacements 
corresponding to N yy and M yy . By use of Hooke's law (2.15) and the 
equilibrium equations (2.1) we find 



(5.9) 


It is convenient to use the following two functions as unknowns. 


9(X) * yVo+ 3x V(X,y) 
fM * & 6 y< Xiy) 


Their integral representations are 


(5.10) 

(5.11) 


g(x) • h 



9 -ix? 

Z Rj(e)> e d? 


(5.12) 


f(x) 


* 


5 X* 




U) + 


-lx? 

5 2 A 1 (5))e d? (5.13) 


The five unknown functions A-j(?) and RjU), j e l,..4 are now elim- 
inated by use of the three homogeneous boundary conditions and inver- 
sion of (5.12) and (5.13). It Is noted that g(x)=f(x)*0 for x^C-l *1 3- 
After a few reductions we find 


, f l it? 

A-jU) a -2 I f(t)e dt 


(5.14) 
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(5.15) 

(5.16) 

(5.17) 

(5.18) 

Due to the systematic nature of these equations the solution Is 
straightforward. When use Is made of the characteristic equation (4.10) 
we get 

■ ^(p 1 -p 2 )(p 1 -p 3 )(pFp 7 

♦ (XX 2 )Y(1-kp,> | g(t)c 1te dt) (5.19) 

The expressions for R 2 (5), R 3 (£) and R^U) are found by interchanging 
the Indices. 

It should be noted that A-|(5)=0 does not follow from the limit 
process k+ 0. Special care must therefore be taken, when relating re- 
sults of the present theory to results from 8th order shallow shell 
theory. 


and the four equations 

I 00,11,(5) • -1 I f(t)e 1t5 dt 
j-1 3 3 il 

* -< 

j, p j'S R j U) c 0 

j, p j'S R j (5) ' 1 1 V)e 1t?dt 
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6. SINGULAR INTEGRAL EQUATIONS 

When A-|(£) and Rj(€) from (5.14) and (5.19) are substituted Into 
(5.2) and (5.4), Integral expressions for N yy (x,y) and M yy (x,y) are 
obtained In terms of the two unknown functions f(t) and g(t). For 
y«0 the Integrals are defined by the limit y+0+. Contributions to the 
Integrals, which arc non-lntegrable for y=0 are extracted by use of 
the equations (5.14)-(5.18). 

The Identity (5.17) Is multiplied by C(?|(A 2 /A) 2 exp(-| 5 y|) and 
added to the Integrand In (5.2). 

NyyU.y) = Effjx* R j (5) 


n+ M e ( V l?|)1 V ,5y, ; iX ' 


j 


d£ 


( 6 . 1 ) 


Asmymptotlc expansion for large values of ? yields 


|F| (nylsDlYl , i Pi 3 Pf . K 

+ m e i - ,. (1 . ...)e 


(m. + l5l)|y| 


( 6 . 2 ) 


Substitution of (6.2) Into (6.1) and use of the identity (5.16) lead 
to a non-lntegrable term of the form 

C 1 f® 1 i ** | ct m 4 “!£yi 

">.o> * $ H ^ V e),e e d5 

i ri -|cy| U(t-x) 

= lim g(t)J sgn(E)e e d£dt 

= > 3 > 
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The last Integral must be evaluated as the Cauchy principal value. 


By a similar, but sightly more complicated procedure, the follow- 
ing non-lntegrable contribution to the moment M yy (x,0) Is found 

MjyU.O) a | ^ dt < 6 - 4 ) 

After extraction of the non-lntegrable parts of the Integrals the 
boundary conditions take the form of two singular Integral equations. 


|^t^x + + j^i 2 (t-x)f(t)dt = 2mN yy (x,0) 


-1<X<1 


(6.5) 


|V 21 (t-x)g{t)dt + [^Tir Fx + K 22 (t-x)]f(t)dt » 2ir — H yy (x,0) 


■l<x<l 


( 6 . 6 ) 


When the symmetry properties m j(5) sm j("5) and PjU)=Pj(-5) are 
noted, the functions ^(t-x) are easily brought In the following form 

w*-*) ■ l> U (p r p 2 )(p r p 3 ' HPi - p ^- i} sin5<t - ,i)d5 l6 - 7) 


r X 2 r 3 Pi + 0*v)C 
K 12 (t-x) = K 21 (t-x) = j o (2( T ') l mi (p r p 2 )( pi -p 3 )(p r p 4 ) } 


s1n?(t-x)d? 


( 6 . 8 ) 



(6.9) 


+ 2k > s1n|{t-x)d? 


The summations imply interchange of the indices. 

!* 

Due to the presence of terms- of order £ 2 the expression (6.9) is 
not suitable for numerical integration. The problem Is solved by con- 
structing the following identity from (5.15), (5.18) and (5.19). 


A 2 "J 2 CPi + 0-v)5 2 ] 2 

Z( T> f pl'(p7-p^TCp7-p 3 Tlp r P 4 ) + 2K 


(1-v) 1 

-TT“ 


£ 2 + 4 



■ 0 


( 6 . 10 ) 


K 22 (t-x) can then be evaluated by the more suitable expression 



+ 2k $U+r 1 ) + > sin5(t-x)d£ (6.11) 


In order to get continuity of the displacements outside the crack 
the solution must also satisfy the two conditions 


J 1 g(t)dt = 0 , |V(t)dt = 0 (6.12) 

The integral equations given above only apply when k>0. In order 
to obtain the corresponding equations for the 8th order theory, where 
k=0, the function A^(i-) must explicitly be set equal to zero. Only 
Myy's dependence on f(t) includes A^f;), and the changes are therefore 
restricted to the last integral in (6.6). The coefficient (l-v 2 )/X H 
must be replaced by (l-v)(3+v)/X 4 , and the kernel K 22 (t-x) is now 
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(6.13) 



[p, + 0-v)S 2 ] 2 

. cp r p 2 Mp r p 3 i(p r p 4 y ' 5+ffl i } 


+ > s1nC(t-x)d? 


7. THE ASYMPTOTIC STRESS FIE LV 


Is 


The solution to the two singular Integral equations (6.5) and (6.6) 
of the form [13] 

g(t) . (l-t 2 rSi(t) 


f(t) . (l-tM'N'tt) 


(7.1) 


In order to obtain the asymptotic stress field In a neighborhood of the 
crack tip X“1 the following formula, which Is derived In Appendix B» Is 
used 



~ {G(l)exp[i(E-sgn(5) J)] 


+ G(-l)exp[-1(£-sgn(C) J)]> 


(7.2) 


Substitution of (7.2) and the similar formula for f(t) Into (5.14)- 
(5.18) makes possible an asymptotic analysis similar to the one which 
lead to the extraction of the Cauchy integrals. Now, however, (6.2) 
must be expressed in the slightly different form 


+ JlL e < ' n j +l5l)|y| (H|5y|)+| OK^HI+Cjlol' 3 ) (7.3) 

j j e 2 
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f 


* 


Cj» j*l,2»... are constants, which do not Influence the result. 

The asymptotic expressions are obtained by application of the 
following formula [14; .3.944] 


c 




fsin] 

cos 


(«5)d5 


jM 


(B 2 +<$ 2 ) 


W 


sin 

cos 


(p arctg |) (7.4) 


It Is seen that a factor (SB)** under the integral sign does not change 
the order of the result for 3»6*0. The order of the terms in (7.3) can 
then be evaluated as if |5y| were a constant. It turns out that only 
the first term of (7.2) contributes to the singular stress field around 
x*l, 

Substitution of (7.2) and (7.3) Into (6.1) yields 

N yy (x,y) ' £ j= (HCMle-^'sInUO-x) -J]d£ <™> 

In the same way we find 

N xx (x.y) - f Q (l-5|y| ) e" 5 l y ls1nU(l-x) - J]dc (7.6) 

N xv (x,y)i~ I—: C 71 ye'^cosUO-x) - J]d? (7.7) 

xy 2/z¥ 

The asymptotic expressions (7. 5)- (7. 7) are also valid for k= 0. For 
. k>0 the asymptotic expressions for the moments are 

M yy(x* y ) ~ ifel^f /f (l+5|y|)e" ?lyl sin[^(l-x)-J]d5 (7.8) 
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M X)l (x.y) - ^ (l-5|y|)e' 5 l y lsin[5(l-x)-J]d5 (7.9) 

V* ,y) ' Wzffi lo ^ ye _e,y| cos[gO-x) -J]d5 (7.10) 

It follows from (7.5)-{7.'i0) that the asymptotic membrane and 
bending stress fields are self similar for k> 0. This is in agreement 

with results reported in [1-4] and [11]. For k=G the asymptotic bend- 

ing stress field is found to be different and depend on Poisson's 
ratio v [7,10]. 

When the coordinates r and 0 are defined by 

x-1 = r cose , y = r sin8 (7.11) 

application of (7.4) and use of trigonometric formulae lead to the 
result 

(7.12) 


For specially orthotropic materials r and 9 are not polar coordinates. 
With the usual definition of the stress intensity factor 

K, = lim /2(X r a) o„{X,,0) (7.13) 

' X-,+3 1 tL 1 

we get 


'xx 


yy 


G(,l)+zF(l )_ 

2/2r 


a 


xy 


^ cos(0/2) + cos (50/2) 
^ cos(0/2) - ^ cos(50/2) 
} sin(0/2) + 1 sin(50/2) 
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(7.14) 


K, * - y E /» [s(l) + zF(l )] 


By Inspection of the Integral equations (6.5) and (6.6) It Is seen that 
the stress Intensity factor at th? surfaces Z*Th/2 depends only on the 
parameters v, X 2 and k* Thus the thickness effect and the Influence of 
the effective transverse shear modulus are combined In the parameter 
K-EABX-). 

: \ 

I. NUMERICAL SOLUTION 


The singular Integral equations (6.5) and (6.6) are solved by use 
of a quadrature formula of closed type developed In [15]. The Integral 
equations are replaced by the algebraic equations 
n 


■ vv> 


( 8 . 1 ) 


^/n.k^l^W^V + t i -x k + K 22^ t 1" x k^ F ^ t 1^ * a" M yy^ x k, 0 ^ 

( 8 . 2 ) 


2h 


where 


- 1-1 


s cos it) i i s l »2 » • - • »n 
x k = cos (f^*) » k=l,2,...,n-l 


(8.3) 

(8.4) 


and 

Vl ■ W n.n * StifclT ,„■! 

(8.5) 

W n,1 “ nh ’ 
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( 8 . 6 ) 


The extra conditions (6.12) are replaced by 

I M .6(0 - 0 . I M .F(t.) - 0 
1«1 n » 1 1 1 >] n » 1 ' 

The convergence of the method Is estimated by evaluation of the 
coefficients bj In the following expansion 

™ 

G(t)«^G(l)+1^6(-l) + (l-t 2 ) l bjUj(t) (8.7) 

where Uj(t) is the Chebyshev polynomial of the second kind of degree j. 
The coefficients bj can be expressed explicitly In terms of the solu- 
tion to the algebraic equations [15], 

b j*n?T CG(t i )-- T -iG(l)-- r iG(-1)]U j (t i ) (8.8) 

The necessary values of the bounded kernels K^tt^-x^), 
and ^(t-i-x^) were calculated by use of Filon's integration formula. 

9. RESULTS 

Numerical results are given in the form of stress Intensity factors 
for two loading situations* constant rumbrane load and constant bending 
moment. For each loading situation two stress Intensity factors are 
given, one for the average stress and one for the surface stress from 
bending. Introduce the following normalization of the stress Intensity 
factors. For compressive membrane load a m 

5 (9.D 
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(9.2) 


k bm‘ CK,(h/2)-K,(0)](a ni /r) 

For a constant bending moment with maximum surface stress o K 


bb - CK,(h/Z) - K,(0)3(a b ^T)" 1 

(9.3) 

Ml * K,«0)(« b 

(9.4) 


Jhe results are given as functions of the shell parameter X 2 for 
various values of a/h in Tables 2-5 and Figs. 3-6. The effective 

c r 

transverse shear modulus yj’p^-has been used together with \f»0.3. 

The first column of the tables and the dashed curve In the figures 
correspond to 8th order shell theory (k*0). As expected from the dif- 
ference in the asymptotic moment fields the bending stresses show some 
differences. The membrane stresses from the 8th order theory, however, 
are found to be representative for very thin shells h/a<10. it is Im- 
portant to note that In general 8th order shell theory gives non- 
conservative estimates of the membrane stresses. 

The extrapolated values of for X 2 =0 are - ' n good agreement 
with the results obtained in [2-4] for a plate. 
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APPENDIX A 


Solution of the equation p* + (X 2 ^)**<1-kp) - 0 
Following the method In [16] the solution of the quartlc equation 
P 4 + (A 2 5) 4 0-kp) * 0 * (A.l) 

requires the real root of the cubic equation 

x 3 - 4(X 2 C)**x - k 2 (X 2 ?) # * 0 (A. 2) 

Now Introduce the parameter 

2 1 

n ■ (-^-) T"‘ (A.3) 

The real solution to (A. 2) Is determined as the sum of Bj and Bg given 
by 

<h’,2 ■ (^)’U 2 5)*[ni^m (A.4) 

i ~) 

This relation Is rewritten in the form 

3 6 (exp (t 1 Arccosn) » 0<n<l 

B l,2 * <*25) «*. " (A.5) 

[exp (+ Arccoshn) » l<n 

The real root of (A. 2) is then given by 

fcos (y Arccosn) > 0<n<l 

X * S,+e 2 • (Xj,5) 2 J 7 . (A-6) 

[cosh Arccoshn) t l<n 

The roots p of (A.l) are found from ! 
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(A. 7) 


p 2 ? li p + ■ 0 

leading to 

* 

P c x 2 5 / co ^a C+ 1(+)1 ^n/iosW5f) 

0<n<l n * cose 


\ 


(A. 8) 


For l<n cos( ) Is replaced by cosh ( ). 


APPENDIX B 

Asymptotic expansion of 1(C) * f 1 4= e 1?t dt 

il /Pt* 

Rewrite the Integral In the form 

HC) 3 j 1 e 15t dt 

-1 

3 ( {[<J»(t)+H-t)]cos(|t)+1[(| ) (t)-<^(-t)]s1n(ct)} -4== (B.l ) 
; o /l-t 2 

(B.2) 


(B.3) 


Tj(t) Is the Chebyshev polynomial of the first kind of degree j» and 
J n (C) is the Bessel function of the first kind of order n. The Integral 
then Is 

oo 

^5) B Tr . c 2j^2j^ I ^! ) + ^ sgn(C)c2j + i d2j + i ( |C| )] (B.4) 

From the asymptotic formula [14] 


Introduce the following series expansion of <f>(t) 

♦<t) 3 I CjTj(t) 
j=0 3 3 

and use the results [14; 7.355] 

j 1 W‘>s1n(5t> ^ = (-lfj0 2n+ ,(|5|) S gn(5) 
0 

£ T 2n(t)cos(«)^.t-l) n |J 2n (|6|) 
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J„(UI) ./^costUI^-^-’n (B.5) 

we get 

\{£) ~ l [c 2;} cos(U| - J) + 1 sgnU)c 2 j +1 sin(|£| -J)] (B.6) 


By use of T 2n (+l)=l and T 2n4l (+!)=+! the leading term of the asymptotic 
expansion of I (£) Is found to be 

1(5) ~ /jfa U0)exp[1(?- sgn(e) J)] 


+ d>(-l )exp[-1 (5 - sgn(t’) J)]} 


(B.7) 


/ 


TABLE 1 

Dimensionless variables. 
Coordinates 


— * £1.1 1 
rt a 

h 



Displacements 



1 u 2 

A * 

-S 


e y = 7T e 2 


j ' 

Stresses and stress resultants 

a - 011 

a x % " ir 

6cT 22 

a yy “ e 

„ . °12 
a xy F“ 

X 

X 

II 

6n 22 
N yy = “RT 

N * Nl2 
N xy "-KT 

M ii 

M x X = 6PF 

H - 6H 22 

”yy “ W 

M - M]Z 
M xy " W 

v = V 1 

* /5hB 

/Sv 2 
v y a ~W~ 

. B F 
* I^FiT 


Parameters 


X 4 = 12(1 -v 2 ) fi2ai * 
■ 1 ' h j R; 

* i2(i-v 2 ) 

x; 2 = 12(1 -v : 

X* = 12(l-v 2 ) 4- 

h 

k = 4tI 




Table 2. Stress intensity factors 


x 2 

tc* 0 

a/h=l 

a/h s 2 

a/h=5 

a/h=10 

0.01 

1.000 

1.000 

1.000 

1.000 

1.000 

0.25 

}.oU 

1.015 

1.015 

1.015 

1.015 

0.50 

1.056 

1.061 

1.058 

1.057 

1.057 

0.75 

1.119 

1.135 

1.123 

1.120 

1.119 


1.198 

1.233 

1.208 

1.200 

1.199 

1.5 

1.391 

1.485 

1.420 

1.398 

1.394 


1.613 

1.788 

1.668 

1.625 

1.618 

UK 

2.095 

2.478 

2.220 

2.122 

2.105 

4.0 

2.588 

3.254 

2.808 

2.634 

2.603 

5.0 

3.075 

4.100 

3.414 

3.146 

3.096 

i m 

3.552 

4.944 

4.069 

3.656 

3.580 

HI 

4.021 

- 

4.723 

4.154 

4.054 

8.0 

4.484 

- 

- 

4.649 

4.515 

10.0 

5.376 

- 

- 

- 

5.422 

12.0 

6.297 

tm 

- 

- 

•* 







Table 3. Stress intensity factors 


^2 

, K =0 

a/h=l 

a/h=2 

j 

a/h=5 

a/h=10 

0.01 





0.0001 

0.25 




1 

0.0208 

0.50 





0.0544 

0.75 

0.142 

0.0951 



0.0891 

1.00 

0.194 

0.1243 


0.1193 

0.1200 

1.5 

mm 

0.1622 

0.1604 

0.1636 

0.1674 

2.0 

nsi 

0.1757 

0.1748 

0.1851 

0.1942 

3.0 

0.371 

0.1507 

0.1397 

0.1661 

0.1887 

4.0 

0.313 

0.0801 

0.0395 

■ 

0.1156 

5.0 

0.176 

-0.0266 

-0.1089 


-0.0140 

6.0 

-0.025 

-0.1510 

-0.2965 


-0.1865 

7.0 

-0.279 

- 


BBB1 


8.0 

-0.579 

- 


-0.7369 

-0.6343 

10.0 

-1.306 

- 


- 

-1.1829 

12.0 

-2.186 

- 


- 

- 
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Table 5. Stress Intensity factors k 


mb 


mi 

o 

n 

a/h=l 

a/h=2 

a/h=5 

a/h=10 

BB 

0.0000 

0.0000 



0.0000 

MM 

0.0043 

0.0069 

1 

0.0051 

0.0047 


0.0115 

0.0184 

0.0158 


0.0126 

0.75 

0.0191 

0.0302 

0.0261 

0.0227 

0.0212 

1.00 

0.0267 

0.0414 

0.0363 


0.0294 

1.5 

0.0407 

0.0607 

0.0544 

0.0478 

0.0447 

WM 

0.0530 

0.0761 

0.0698 

0.0619 

0.0582 

■9 

0.0731 

0.0977 

0.0933 

0.0846 

0.0799 

111 

0.0882 

0.1121 

0.1096 

0.1012 

0.0961 

5.0 

0.0998 

0.1223 

0.1209 

0.1136 

0.1081 

6.0 

0.1088 

0.1282 

0.1303 

0.1231 

0.1175 

7.0 

0.1161 

- 

0.1370 

0.1303 

0.1249 

8.0 

0.1221 

- 

- 

0.1359 

0.1308 

10.0 

0.1309 

- 

- 

- 

0.1399 

12.0 

0.1388 

- 

- 

- 

- 
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Fig. 5. 


Stress intensity factor 



